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Electron-impact  dissociative  ionization  of  ethylene 
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The  ionization  rates  of  the  electron-impact  ionization  and  dissociative  ionization  of  ethylene  for  two  typical 
cases:  thermal  plasmas  with  an  electron  temperature  in  the  range  between  1000  and  24  000  K;  and  nonthermal 
plasmas  with  the  reduced  electric  field  in  the  range  between  10  and  200  Td  are  presented.  Electron-impact 
dissociative  ionization  rates  were  calculated  for  11  fragment  ions  of  ethylene  in  the  case  of  low  temperature 
thermal  plasma,  and  in  a  case  of  nonthermal  ionized  mixture  of  argon  and  ethylene.  Dissociative  ionization 
cross  sections  were  calculated  using  a  semiempirical  binary-encounter  bethe  (BEB)  model  [Y.  K.  Kim  and  M. 
E.  Rudd,  Phys.  Rev.  A  50,  3954  (1994)],  with  each  of  the  four  most  dominant  fragments,  C2H4+,  C2H3+, 
C2H2+,  and  H+,  being  associated  with  a  single  molecular  orbital.  Calculated  cross  sections  are  used  in  this 
analysis  due  to  the  fact  that  the  existing  experimental  data  are  the  least  accurate  in  the  threshold  region  and  the 
calculated  results  can  improve  the  accuracy  in  that  region.  Also,  the  procedure  may  be  extended  to  molecules 
for  which  experimental  data  are  not  available. 

DOI:  10. 1 103/PhysRevA.73.0227 1 1  PACS  number(s):  34.80.Gs,  34.80.Ht,  52.20.Fs,  52.20.Hv 


I.  INTRODUCTION 

Chemical  kinetic  models  of  ignition  and  combustion  un¬ 
der  plasma  conditions  require  reliable  data  on  total  and  on 
dissociative  ionization  processes.  Total  electron-impact  ion¬ 
ization  rates  are  known  rather  accurately  for  atmospheric  or 
small  hydrocarbon  molecules,  but  very  few  data  exist  on 
more  complex  hydrocarbons  or  on  dissociative  ionization 
rates  and  the  production  of  hydrocarbon  fragment  ions.  In 
the  case  of  more  complex  hydrocarbons,  where  the  experi¬ 
mental  database  is  limited  to  a  few  molecular  species,  calcu¬ 
lation  techniques  with  advanced  predictability,  taking  into 
account  the  current  knowledge  of  the  collision  cross  sections, 
are  very  much  needed.  In  addition,  ionized  fragments  of  eth¬ 
ylene  are  found  in  the  interstellar  space  and  in  the  planetary 
ionospheres.  However,  our  goal  is  to  understand  the  disso¬ 
ciative  ionization  process  in  low  electron  energy  regimes  of 
electric  discharge  plasmas. 

Dissociative  ionization  (DI)  processes  are  very  important 
in  plasmas  and  DI  processes  involving  hydrocarbons  will  be 
very  important  for  new  applications  of  plasma  to  enhance 
combustion.  Rates  for  the  generation  of  radical  ions  by  the 
electron  impact  on  hydrocarbon  molecules  are  an  important 
subset  of  data  for  the  chemical  kinetic  modeling  of  a  wide 
range  of  the  plasma-assisted  processes  involving  hydrocar¬ 
bons.  Besides  the  obvious  impact  on  ion  and  neutral  compo¬ 
sition,  radical  ion  generation  processes  define  charge  particle 
balance  and  electron  temperature  in  molecular  plasmas. 
These  data  are  largely  unavailable  and  are  usually  derived 
indirectly,  by  elimination,  rather  than  by  experiment  or  cal¬ 
culation.  One  of  the  problems  with  the  current  experimental 
data  on  fragmentary  ionization  processes  is  that  the  practical 
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application  usually  requires  a  behavior  of  cross-section  data 
at  low  energies,  near  the  threshold,  where  some  of  the  avail¬ 
able  experimental  methods  are  the  least  accurate  [1].  More 
accurate  experimental  techniques  [2]  have  not  been  used  to 
study  ethylene.  On  the  other  hand,  recently  developed  calcu¬ 
lations  give  more  or  less  correct  values  of  the  total  ionization 
cross  section  [3-5],  but  are  still  unable  to  produce  full  data 
sets  for  fragment  ion  generations. 

In  this  paper,  the  quantitative  analysis  of  the  electron- 
impact  ionization  of  ethylene,  where  the  partial  ionization 
cross  sections  for  fragment  generation  are  correlated  with  the 
molecular  structure  is  presented.  Successful  completion  of 
this  analysis  may  lead  to  new  semiempirical  calculation 
models  that  combine  the  experimental  data  on  molecular 
structure  and  fragmentation  processes  with  the  binary- 
encounter-bethe  model  [3-5]  to  produce  reliable  fragmentary 
ionization  rates  for  an  application  in  discharges  involving 
hydrocarbons. 


II.  TOTAL  CROSS  SECTION  FOR  ELECTRON-IMPACT 
IONIZATION  OF  ETHYLENE 


Kim  and  co-workers  [3-7]  have  provided  a  useful  formula 
for  electron-impact  ionization  using  a  theory  referred  to  by 
these  authors  as  the  binary-encounter-bethe  (BEB)  model. 

The  basic  formula  for  the  ionization  cross  section  per  mo¬ 
lecular  orbital  is  [7], 


&BEB  -  ' 


t  + 


u  +  1 


+  (2 -Q) 


In  t  \ 
t  +  1  / 


(i) 

where  t=T/B,  u=U/B,  S=4pa02N(R/B)2,  a0= 0.529  18  A, 
R  =13.6057  eV,  B  is  the  binding  energy,  U  is  the  orbital  ki¬ 
netic  energy,  N  is  the  electron  occupation  number,  Q  is  the 
dipole  constant,  n- 1  for  neutral  molecules,  and  n-  2  for  the 
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TABLE  I.  HOMO  binding  energy  and  ionization  energy  for 
C4H4+.  Experimental  value  for  the  ionization  energy  of  C2H4  is 
10.51  eV. 


Model 

HOMO  binding 
energy  (eV) 

Ionization 
energy  (eV) 

RHF/6-31G(d) 

10.19 

8.78 

RHF/6-311G(2rf/) 

10.37 

8.90 

B3LYP/6-3 1 1 +G(3  df,  2 p) 

7.69 

10.33 

G1 

10.21 

10.59 

G2 

10.26 

10.58 

G2(MP2) 

10.26 

10.59 

singly  charged  molecular  ion.  The  dipole  constant  Q  is  de¬ 
fined  [6]  in  terms  of  the  continuum  dipole  oscillator  strength 
df/dW , 

2  f  B  df 

Q=-\ - —dW,  (2) 

NJ  B  +  WdW 

where  W  is  the  kinetic  energy  of  the  ionized  electron.  By 
assuming  that  a  simple  form  for  df/dW ,  namely,  df/dW 
which  approximately  describes  the  asymptotic 
behavior  of  the  oscillator  strength  in  hydrogen  atoms  in  the 
ground  state,  reliable  differential  ionization  cross  sections 
can  be  obtained  for  molecules  with  simple  shell  structures 
including  highly  symmetrical  hydrocarbon  molecules,  such 
as  ethylene.  In  addition,  this  assumption  results  in  Q=  1 
which  greatly  simplifies  the  calculation  of  the  cross  section. 

The  total  electron-impact  ionization  cross  section  is  then 
obtained  by  summation  over  all  occupied  molecular  orbitals. 
In  order  to  obtain  the  parameters  relating  to  the  molecular 
orbitals,  ab  initio  quantum  mechanical  calculations  are  per¬ 
formed.  One  difficulty  encountered  in  performing  these  cal¬ 
culations  for  hydrocarbon  ions  is  the  accurate  calculation  of 
low  lying  electronic  states;  however,  the  BEB  model  only 
requires  information  regarding  the  ground  state  wave  func¬ 
tion.  Calculations  have  been  carried  out  at  the  Hartree-Fock 
level  of  theory  with  Gaussian  basis  sets  [8]  in  order  first  to 
calculate  the  molecular  geometry.  The  resulting  wave  func¬ 
tions  were  varied  to  determine  if  it  corresponded  to  the 
lowest-lying  electronic  state — and  if  not,  the  multiplicity  or 
electron  configuration  was  adjusted  so  as  to  yield  the  ground 
state. 

Higher-level  calculations  including  the  G2(MP2)  method 
of  GAUSSIAN  98  [8]  were  also  performed  but  did  not  change 
the  results  appreciably  (see  Table  I)  implying  that  the  use  of 
the  higher-level  calculations  are  not  warranted  in  the  case  of 
the  highest  occupied  molecular  orbital  (HOMO)  binding  en¬ 
ergy.  The  observation  that  Hartree-Fock  wave  functions  are 
adequate  has  been  observed  for  many  molecules  [3-5,7].  The 
orbital  parameters  B ,  U,  and  N  are  obtained  by  the  quantum 
chemistry  calculations  and  are  listed  in  Table  II. 

The  molecular  orbital  constants  N ,  B ,  and  U  can  easily  be 
calculated  using  standard  molecular  wave  function  codes. 
However,  BEB  cross  sections  at  low  temperature  are  sensi¬ 
tive  to  the  lowest  binding  energy  and  the  experimental  ver- 


TABLE  II.  Molecular  orbitals  (MO),  electron  binding  energy 
(5),  orbital  kinetic  energy  (£/),  and  electron  occupation  number  (V). 
B  and  U  values  are  obtained  using  the  G2  method. 


MO 

B  (eV) 

U  (eV) 

N 

i«g 

305.64 

435.74 

2 

l*lu 

305.59 

436.19 

2 

2  ag 

28.12 

40.52 

2 

2bla 

21.59 

33.55 

2 

1*2„ 

17.48 

25.23 

2 

3ag 

16.03 

34.71 

2 

l*3g 

13.80 

28.59 

2 

1*3» 

10.26 

25.83 

2 

deal  ionization  potential  should  be  used  if  available,  unless 
the  calculated  value  of  the  binding  energy  has  included  the 
effect  of  electron  correlation.  Therefore,  the  average  mea¬ 
sured  ionization  potential  for  ethylene  and  the  appearance 
potentials  for  fragments  ions  have  been  presented  in  Table 
III.  The  mean  fraction  of  the  total  cross  section  was  taken 
over  the  energy  range  10-200  eV. 

Total  ionization  cross-section  data  for  ethylene  are  given 
in  Fig.  1  for  the  energy  range  that  is  relevant  for  plasma 
kinetic  studies,  that  is  from  the  threshold  energy  to  100  eV. 
The  two  sets  of  results  calculated  using  the  HOMO  binding 
energy  and  using  the  experimental  ionization  potential,  are  in 
fairly  good  agreement.  The  discrepancy  between  calculated 
and  experimental  data  is  still  substantial.  Early  measure¬ 
ments  of  absolute  total  ionization  cross  section  by  Rapp  and 
Englander-Golden  [10]  are  in  good  agreement  with  calcula¬ 
tions  around  the  threshold  energies,  but  higher  by  about  25% 
at  higher  energies  (>40  eV).  Experimental  results  of  Nish- 
imura  and  Tawara  [11]  are  in  rather  good  agreement  with 
calculations  over  the  full  energy  range,  especially  near  the 


TABLE  III.  Mean  fraction  of  the  total  cross  section  for  the  ob¬ 
servable  fragment  ions  of  C2H4.  (*  obtained  by  electron-impact 
method  exclusively). 


Fragment 

ion 

Appearance 
potential  (eV) 
[6]*,  [9] 

Evaluated 

appearance 

potential 

(eV) 

Mean 

fraction  of 

total  cross- 
section  [1] 

Relative 

uncertainty 

c2h4+ 

10.51 

10.51 

0.336 

±6% 

c2h3+ 

13.10-13.52 

13.06 

0.233 

±9% 

C2H2+ 

13.14 

13.20 

0.212 

±3% 

c2h+ 

18.70-19.20 

19.22 

0.039 

±24% 

c  + 

l_2 

24.50 

24.06 

0.013 

±26% 

ch4+ 

18.66 

18.82 

ch3+ 

16.95 

16.96 

0.0025 

±17% 

ch2+ 

17.82 

17.95 

0.0482 

±13% 

CH+ 

17.68 

17.76 

0.0201 

±31% 

c+ 

18.30 

17.37 

0.0129 

±28% 

H2+ 

22.40 

23.04 

0.0054 

±26% 

H+ 

18.66 

18.41 

0.0773 

±28% 
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Electron  Energy  (eV) 


FIG.  1.  Total  cross  section  for  electron-impact 
ionization  of  ethylene.  Experiment:  ♦,  Ref.  [1]; 
□,  Ref.  [11];  A,  Ref.  [10];  •,  Ref.  [12].  Calcu¬ 
lation:  thick  solid  line,  BEB  calculation  with  the 
binding  energies  of  the  highest  three  OMO  sub¬ 
stituted  by  the  appearance  potentials  of  C2H4+, 
C2H3+,  and  C2H2+,  respectively;  thin  solid  line, 
calculation  using  computed  binding  energies  and 
experimental  ionization  energy  for  HOMO. 


threshold.  The  discrepancy  is  reduced  by  substituting  the 
binding  energy  of  several  higher  occupied  molecular  orbitals 
with  the  appropriate  experimental  appearance  potentials  for 
certain  ion  fragments.  Using  this  approach,  the  agreement 
between  the  experimental  and  calculated  data  below  70  eV  is 
achieved  with  the  overall  uncertainty  of  5%. 

Values  for  the  absolute  total  cross  section  at  high  energies 
obtained  by  Tian  and  Vidal  [1]  are  in  a  good  agreement  with 
other  data  in  the  energy  range  between  25  and  60  eV.  The 
values  near  the  threshold  are  substantially  lower  due  to  a 
experimental  normalization  procedure  based  on  the  argon  ion 
counting.  Ethylene  fragment  ions  with  a  low  appearance  po¬ 
tential  (binding  energy)  were  not  possible  to  observe  in  that 
experiment  because  of  the  high  value  of  argon  ionization 
potential,  which  is  15.756  eV.  Except  for  the  systematic  er¬ 
ror  near  the  threshold  for  ions  with  lower  appearance  poten¬ 
tials,  data  at  energies  around  50  eV  give  very  good  founda¬ 
tion  for  the  analysis  of  ethylene  ion  fragmentation.  Finally, 
the  value  of  5.86  X  10-16  cm2  at  70  eV,  obtained  using  an 
additive  formula  of  Ref.  [12]  is  also  in  a  fairly  good  agree¬ 
ment  with  the  experiment  and  calculation. 

III.  PARTIAL  IONIZATION  CROSS  SECTIONS 

In  this  section  the  possibility  to  predict  the  fragmentation 
branching  fractions  based  on  the  molecular  structure  is  ex¬ 
plored.  This  section  outlines  a  process  to  reproduce  the  ex¬ 
perimental  results  on  fragmentation  by  molecular  structure 
reactant-product  analysis.  If  that  exercise  proves  successful, 
we  could  proceed  with  fragmentation  analysis  for  heavier 
hydrocarbons. 

Tian  and  Vidal  [1]  have  published  tabulated  values  of  dis¬ 
sociative  ionization  cross  sections  for  the  electron  impact  on 
ethylene  in  the  energy  range  between  17.5  and  600  eV. 
There  are  three  distinctive  groups  of  12  possible  ionic  frag¬ 
ments,  C2H+  (n  =  0,  1,  2,  3,  and  4),  CH+  (n  =  0,  1,  2,  3,  and  4), 
and  H*  (n=  1,2).  Tian  and  Vidal  [1]  were  able  to  observe  11 
of  them.  In  Table  III  are  given  their  appearance  potentials 
compiled  from  ion  energetics  reference  data  [6],  and  relative 
contributions  of  the  particular  ion  products  evaluated  from 
the  experimental  data  [1]  over  the  electron  energy  range 
50-600  eV. 


An  inspection  of  Table  III  reveals  that  one-third  of  the 
total  cross  section  belongs  to  the  nondissociative  ionization 
by  electron  impact 

£  +  C2H4  — »  C2H4+  +  £  +  £  .  (3) 

Most  probable  experimental  electron-impact  dissociative 
ionization  processes  produce  fragment  ions  without  breaking 


carbon  bonds,  namely, 

£  +  C2H4  — »  C2H3+  +  H  +  £  +  £ 

(4) 

and 

c  +  C2H4  — »  C2H2+  +  H2  +  £  +  £  . 

(5) 

Taken  together,  the  electron-impact  ionization  processes 
leading  to  the  group  C2H+  account  for  more  than  80%  of  ion 
products.  Other  important  processes  are 

£  +  C2H4  — >  H+  +  C2H3  +  £  +  £ 

(6) 

and 

£  +  C2H4  — >  CH2  +  CH2  +  £  +  £  . 

(7) 

Only  reaction  (7)  involves  carbon  bond  breaking.  It  is  inter¬ 
esting  to  note  that  fragment  ions  CH4+  and  CH3+  have  very 
low  ionization  cross  sections,  since  their  formation  requires 
carbon  bond  breaking  and  redistribution  of  hydrogen  atoms. 
CH4+  is  not  observed  in  the  experiment  [1],  in  spite  of  its 
relatively  low  appearance  potential.  We  can  conclude  that 
dominant  fragment  ions  are  those  that  do  not  require  carbon 
bond  breaking. 

A.  C2Hn+  ions 

Since  the  formation  of  the  ethylene  ion,  C2H4+,  does  not 
involve  any  bond  breaking,  it  is  natural  to  assume  that  the 
removed  electron  was  associated  with  HOMO.  Therefore,  we 
first  inspect  the  relation  between  the  experimental  nondisso¬ 
ciative  ionization  cross  section  for  electron  impact  on  C2H4 
and  the  calculated  partial  ionization  cross  section  associated 
with  the  HOMO  with  its  binding  energy  substituted  by  the 
vertical  ionization  potential  IP=  10.51  eV.  The  results  are 
given  in  Fig.  2  for  the  electron  energy  from  the  threshold  to 
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Electron  Energy  (eV) 


FIG.  2.  Partial  ionization  cross  section  for 
C2H4+  formation.  Experiment:  ♦,  Ref.  [1].  Cal¬ 
culation:  thick  solid  line,  calculation  based  only 
on  HOMO  with  the  binding  energy  equal  to  the 
vertical  ionization  potential  (AP=  10.51  eV);  thin 
solid  line,  calculation  using  computed  binding 
energy  for  HOMO. 


100  eV.  The  agreement  is  better  than  5%  at  energies  above 
50  eV,  which  is  comparable  to  the  uncertainties  of  calculated 
values,  mostly  due  to  the  uncertainty  of  the  binding  energy 
(2%)  and  uncertainty  of  the  orbital  kinetic  energy  U  (2.5%), 
based  on  more  than  a  dozen  models  tested  up  to  the  level  of 
G2.  At  energies  from  the  threshold  to  50  eV  we  have  more 
reason  to  believe  that  HOMO  ionization  cross-section  values 
are  more  reliable  than  the  experimental  data,  since  the  ex¬ 
periment  cannot  record  the  threshold  values,  and  the  experi¬ 
mental  error  remains  high  up  to  about  50  eV.  It  has  been 
mentioned  earlier  that  the  experiment  [1]  could  not  repro¬ 
duce  the  threshold  values  correctly,  since  the  appearance  po¬ 
tential  for  C2H4+  is  lower  than  the  ionization  potential  of 
argon,  which  served  as  the  carrier  gas. 

Results  for  other  two  important  ion  fragments  from  the 
same  group,  C2H3+,  and  C2H2+,  are  given  in  Figs.  3  and  4, 
respectively.  Figure  3  shows  results  of  the  calculations  of 
partial  ionization  cross  section  for  the  1  b3g  orbital  together 
with  experimental  data  for  the  fractional  ionization  cross  sec¬ 
tion  for  the  C2H3+  ion  fragment.  The  cross  sections  calcu¬ 
lated  using  the  orbital  parameters  obtained  at  the  HF/6- 
31G(d)  level  of  theory  yielding  for  the  1  b3g  MO  a  binding 
energy  B  of  13.66  eV  are  in  good  agreement  with  the  experi¬ 
mental  data  at  an  energy  above  50  eV.  Higher  levels  of 
Hartree-Fock  techniques  give  slightly  higher  values  of  bind¬ 
ing  energy,  up  to  13.80  eV.  When  the  binding  energy  for  this 
orbital  is  substituted  by  the  average  experimental  appearance 
potential,  AP=  13.33  eV  for  C2H3+,  the  corresponding  partial 
ionization  cross  section  is  slightly  increased.  Taking  into  ac¬ 


count  the  uncertainties  in  the  experimental  and  calculated 
data,  the  two  calculations  can  be  considered  almost  as  iden¬ 
tical.  The  agreement  between  the  two  calculations  suggests 
that  the  partial  ionization  cross  section  for  the  1  b3g  orbital  is 
an  adequate  approximation  for  the  dissociative  ionization 
leading  to  vinyl  ion  C2H3+. 

Quantitative  agreement  between  the  empirical  approach 
and  MO  calculation,  starts  to  break  down  at  the  next  occu¬ 
pied  orbital,  3 ag.  Remedy  for  this  particular  case  is  to  sub¬ 
stitute  the  binding  energy  of  the  orbital  3 ag,  with  the  appear¬ 
ance  potential  for  the  C2H2+,  as  in  the  many  cases  of  total 
ionization  cross  sections  [6].  As  seen  in  Fig.  4,  this  approach 
gives  a  quite  accurate  reproduction  of  cross  sections  both  in 
the  threshold  region  and  at  higher  energies,  where  it  could  be 
tested  by  the  experimental  data.  Partial  ionization  cross  sec¬ 
tions  corresponding  to  the  next  available  molecular  orbital, 
3 ag,  is  about  40%  lower  than  the  experimental  dissociative 
ionization  cross  section  of  the  C2H2+  ion  fragment.  The  large 
difference  between  the  calculated  binding  energy  of  this  or¬ 
bital  (15.97  eV)  and  the  average  experimental  appearance 
potential  (13.17  eV)  is  responsible  for  the  large  differences 
shown  in  Fig.  4  between  the  two  calculations. 

B.  Hydrogen  ions 

A  similar  situation  is  found  in  the  case  of  a  H+  fragment 
ion,  as  seen  in  Fig.  5,  where  the  discrepancy  between  em¬ 
pirical  calculation  and  the  experiment  is  more  profound, 


FIG.  3.  Partial  ionization  cross  section  for  the 
production  of  C2H3+  fragment  ion.  Experiment: 
♦  ,  Ref.  [1].  Calculation:  thick  solid  line,  BEB 
based  on  the  same  MO  with  AP=  13.06  eV  in¬ 
stead  of  binding  energy  of  13.80  eV;  thin  solid 
line,  BEB  based  only  on  the  MO  1  b3g. 
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FIG.  4.  Partial  ionization  cross  section  for  the 
production  of  C2H2+  fragment  ion.  Experiment: 
♦  ,  Ref.  [1].  Calculation:  thick  solid  line,  BEB 
based  only  on  the  MO  3  a  lg  with  an  appearance 
potential  of  13.20  eV;  thin  solid  line,  BEB  based 
on  MO  3a lg  with  a  calculated  binding  energy  of 
16.03  eV. 


but  still  rather  close  to  the  frame  of  experimental  error.  How¬ 
ever,  the  use  of  the  binding  energy  for  the  1  b2n  MO  from 
Table  II  would  lead  to  a  large  overestimate  of  this  cross 
section.  Therefore,  the  binding  energy  for  the  same 
hydrogen-bonding  orbit  1  b2u  was  substituted  by  the  appear¬ 
ance  potential  of  H+. 

C.  Other  ions 

Cross  sections  for  dominant  ion  species  C2H4+,  C2H3+, 
C2H2+,  and  H+  can  each  be  associated  with  a  single  molecu¬ 
lar  orbital  and  calculated  using  the  values  of  appearance 
potentials  instead  of  the  binding  energy.  These  cross  sections 
are  then  in  very  good  agreement  with  the  experiment.  How¬ 
ever,  this  approach  is  not  sufficient,  since  it  is  obvious 
from  Tables  II  and  III  that  there  cannot  be  a  one-on-one 
correspondence  between  the  molecular  orbitals  and  ion 
fragments.  Moreover,  calculated  binding  energies  differ 
considerably  from  the  appearance  potentials  of  particular  ion 
fragments. 


The  sum  of  dissociative  ionization  cross  sections  for  four 
main  fragment  ions,  C2H4+  (see  Fig.  2),  C2H3+  (see  Fig.  3), 
C2H2+  (Fig.  4),  and  H+  (Fig.  5),  calculated  using  their  ap¬ 
pearance  potentials  instead  of  the  calculated  binding  energies 
of  the  assigned  MO,  constitute  more  then  80%  of  the  total 
ionization  cross  section  for  ethylene.  The  remaining  20%  is 
distributed  between  all  other  fragment  ions.  In  Fig.  6,  we 
compare  the  sum  of  BEB  cross  sections  of  the  next  two 
occupied  orbits  of  C2H4,  2MU,  and  2 ag  with  the  sum  of  the 
experimental  cross  sections  of  all  other  observable  ions, 
C2H+,  C2+,  CH3+,  CH2+,  CH+,  C+,  and  H2+.  Here  the  ethynyl 
radical  ion  C2H+,  dicarbon  ion  C2+,  and  hydrogen  molecular 
ion  were  associated  with  the  2 biu  orbital,  and  the  methyl 
radical  ion  CH3+,  methylene  ion  CH2+,  methylidyne  ion 
CH+,  and  carbon  ion  C+,  were  associated  with  the  2 ag  or¬ 
bital.  Appearance  potentials  for  these  ion  fragments  are 
higher  than  the  ionization  potential  of  argon,  and  the  experi¬ 
mental  data  from  Ref.  [1]  should  be  accurate  within  the 
stated  error  of  10%.  This  is  reflected  in  the  fact  that  the 
present  calculation  and  experiment  [1]  are  in  much  better 
agreement  near  the  threshold. 


FIG.  5.  Partial  ionization  cross  section  for  H+ 
fragment  ion.  Experiment:  ♦.  Ref.  [1].  Calcula¬ 
tion:  thick  solid  line,  BEB  using  only  MO  2 ag 
with  AP=22.3  eV  instead  of  binding  energy;  thin 
solid  line,  BEB  based  on  MO  2 ag  with  a  calcu¬ 
lated  binding  energy  of  28.12  eV. 


Electron  Energy  (eV) 
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FIG.  6.  Total  ionization  cross  section  of  the 
remaining  observable  fragment  ions  (CH3+, 
CH2+,  CH+,  C+,  C2H+,  C2+,  and  H2+).  Experi¬ 
ment:  ♦  ,  Ref.  [1].  Calculation:  thin  solid  line, 
BEB  using  the  2 bXxx  and  2 ag  MO  with  corre¬ 
sponding  appearance  potentials. 


IV.  IONIZATION  RATES 

While  the  near- threshold  value  cross  sections  are  not  par¬ 
ticularly  meaningful  in  the  astrophysical  plasmas  where  the 
ionization  processes  are  carried  on  by  more  energetic  par¬ 
ticles,  they  are  very  important  in  the  laboratory  processing 
plasmas.  Since  the  electrons  capable  to  ionize  are  in  the  tail 
of  their  energy  distribution,  the  threshold  energy  and  abso¬ 
lute  value  of  the  cross  section  for  each  fragment  ion  are 
critical  for  determining  the  ionization  rates  in  typical  plas¬ 
mas.  In  evaluating  the  ionization  rates  for  the  dissociative 
ionization  of  ethylene  in  discharge  plasmas,  two  distinct  situ¬ 
ations  are  considered.  The  first  case  assumes  that  the  elec¬ 
trons  in  the  discharge  are  in  equilibrium,  with  their  tempera¬ 
ture  defined  by  the  Maxwellian  energy  distribution.  The 
second  case  assumes  the  discharge  in  a  lean  mixture  of  eth¬ 
ylene  and  argon,  where  electron  energy  distribution  function 
is  defined  by  the  reduced  electric  field  E/N ,  in  the  discharge 
dominated  by  argon  collision  processes.  Both  cases  are  char¬ 
acterized  by  relatively  low  average  electron  energy,  order  of 
1  to  2  eV,  and  are  therefore  sensitive  to  the  near-threshold 
values  of  ionization  cross  sections.  Both  cases  are  typical  of 
laboratory- scale  plasma. 

Ionization  rates  for  thermal  plasmas  are  given  in  Table  IV 
and  have  been  obtained  by  convoluting  the  computational 
cross  section  data  with  a  Maxwellian  electron  distribution 
fuction.  For  the  electron  temperature  range  from 
1000  to  24  000  K,  the  ionization  rate  coefficients  are  well 
described  by  an  Arrhenius  form. 

In  the  more  frequent  case  of  nonthermal  plasmas,  the 
electron  temperature  has  a  somewhat  ambiguous  meaning  of 
the  average  electron  energy,  and  a  more  appropriate  param¬ 
eter  defining  a  particular  electron  energy  distribution  func¬ 
tion  is  the  reduced  electric  field  E/N.  Ionization  rates 
of  fragment  ions  at  10,  50,  100,  and  200  Td  (1  Td 
=  10“ 17  V  cm2)  in  a  lean  mixture  of  ethylene  with  argon,  are 
given  in  Table  V.  It  is  interesting  to  note  that  the  dominance 
of  the  fragment  ions  with  carbon  bonds  preserved  is  much 
more  pronounced.  All  other  fragment  ions  appear  in  traces  at 
low  E/N.  At  high  E/N,  the  appearance  of  C2H+,  CH2+,  and 
H+  can  be  counted  in  parts  of  a  percent. 

The  experimental  cross  section  values  of  the  most  abun¬ 
dant  fragment  ions  [1]  near  threshold  cannot  be  used  due  to 


the  systematic  error  for  product  channels  with  threshold  en¬ 
ergies  below  the  ionization  potential  of  argon  as  previously 
discussed.  Therefore,  in  order  to  compare  the  ionization  rates 
obtained  from  the  BEB  model  calculation  to  those  derived 
from  the  experimental  data,  the  experimental  points  for  elec¬ 
tron  energies  below  25  eV  were  eliminated  and  approxi¬ 
mated  by  linearly  interpolating  the  value  of  the  cross  section 
from  25  to  10.51  eV.  The  maximum  difference  between  the 
ionization  rate  coefficients  derived  from  the  modified  experi¬ 
mental  and  calculated  cross  sections  is  about  35%  at  high 
temperatures  (above  10  000  K)  and  high  E/N  (above  10  Td). 
Based  on  these  results,  the  values  of  ionization  recorded  in 
Tables  IV  and  V  are  accurate  to  35%  for  the  ions  with  ap¬ 
pearance  potentials  below  the  ionization  potential  of  argon, 
i.e.,  the  major  product  channels  C2H4+,  C2H3+,  and  C2H2+. 

For  fragment  ions  with  appearance  potentials  above  the 
ionization  potential  of  argon,  the  experimental  values  near 
threshold  are  reliable.  In  these  cases,  the  ionization  rates 
were  calculated  from  the  experimental  data  using  a  recurrent 
polynomial  formula  in  the  threshold  region  to  make  the  data 
continuous.  This  procedure  results  in  differences  between  the 
ionization  rate  coefficients  derived  from  the  modified  experi- 


TABLE  IV.  Fitted  Arrhenius  parameters  ( k=AeBITe )  for  partial 
ionization  rates  of  ethylene  fragment  ions  in  thermal  discharge 
plasmas. 


Fragment  ions 

A 

(10-8  cm3  s-1) 

B 

(1000  K) 

c2h4+ 

1.7107 

125.47 

c2h3+ 

1.0514 

158.12 

c2h2+ 

0.9632 

156.18 

c2h+ 

0.1661 

225.01 

c2+ 

0.02716 

308.88 

ch3+ 

0.009342 

204.36 

ch2+ 

0.1771 

233.10 

CH+ 

0.08087 

264.67 

c+ 

0.009213 

214.91 

H2+ 

0.01985 

262.80 

H+ 

0.3957 

260.53 
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TABLE  V.  Partial  ionization  rates  of  ethylene  fragment  ions  in  a  very  lean  mixture  of  ethylene  with 
weakly  ionized  argon  in  units  of  (cm1 2 3 4 5 6 7 8  s_1)  for  different  E/N  value. 


Fragment  ions 

10  (Td) 

50  (Td) 

100  (Td) 

200  (Td) 

c2h4+ 

9.63  X  10-14 

8.02  X  10-13 

1.62  X10-12 

4.78  X  10-12 

c2h3+ 

9.04  X  10-17 

4.51  X  10-14 

2.27  X  10-13 

1.09  X  10-12 

c2h2+ 

1.36  X  10-16 

4.82  X  10-14 

2.28  X  10“13 

1.06  X  10-12 

c2h+ 

2.11  X  10-30 

1.14X  10-17 

9.68  X  10“16 

1.75  X  10~14 

c2+ 

8.93  X  10-23 

1.44  X10-18 

2.53  X  10“16 

ch3+ 

5.42  X10-28 

3.35  X  10-18 

1.44  X10-16 

1.90  X  10~15 

ch2+ 

2.68  X  10-27 

3.69  X10-17 

1.73  X10-15 

2.39  X  10~14 

CH+ 

1.86  x  10-19 

9.59  X10“17 

3.44  X  10~15 

c+ 

1.04  X  10-30 

1.77  X  10-18 

1.07  X10-16 

1.82  X  10~15 

h2+ 

1.50  X  lO-20 

1.11  X10-17 

4.84  X  10“16 

H+ 

3.77  X10-19 

2.59  X  IQ’16 

1.16X  10-14 

mental  and  calculated  cross  sections  is  about  50-70  %  for 
the  ionization  rates  for  the  minor  product  channels. 

V.  CONCLUSION 

Ionization  rates  of  electron-impact  ionization  and  disso¬ 
ciative  ionization  of  ethylene  for  two  typical  cases:  thermal 
plasmas  with  electron  temperature  in  the  range  between  1000 
and  24  000  K;  and  nonthermal  plasmas  with  the  reduced 
electric  field  in  the  range  between  10  and  200  Td  have  been 
calculated.  The  values  obtained  for  ionization  rates  are  based 
on  the  calculated  absolute  dissociative  ionization  cross  sec¬ 
tions  using  the  binary-encounter-bethe  model  where  single 
MO  for  dominant  fragment  ions  was  assumed  and  assigned 
but  the  corresponding  binding  energy  was  substituted  by  the 


appearance  potential  for  the  particular  ion  fragment.  The  cal¬ 
culated  values  have  been  compared  favorably  with  available 
experimental  results.  The  sum  of  the  fragment  dissociation 
cross  section  improves  the  agreement  between  calculated 
values  and  experiment.  Presented  results  demonstrate  the  in¬ 
crease  of  dominance  of  fragment  ions  with  carbon  bonds 
preserved  in  low-temperature  laboratory  plasma,  due  to  their 
lower  ionization  threshold  and  higher  absolute  cross  sec¬ 
tions.  Present  empirical  approach  to  the  calculation  is  going 
to  be  further  tested  for  more  complex  hydrocarbons. 
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